Electronic Inhomogeneity and Breakdown of the Universal Thermal Conductivity in 

Cuprate Superconductors 
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We report systematic, high-precision measurements of the low-T (down to 70 mK) thermal con- 
ductivity K of YBa2Cu30y, La2-iSra;Cu04 and Bi2Sr2CaCu2 084-i • Careful examinations of the 
Zn- and hole-doping dependences of the residual thermal conductivity kq/T, as well as the in-plane 
anisotropy of ko/T in Bi2Sr2CaCu208+i , indicate a breakdown of the universal thermal conductiv- 
ity, a notable theoretical prediction for d-wave superconductors. Our results point to an important 
role of electronic inhomogeneities, that are not considered in the standard perturbation theory for 
thermal conductivity, in the under- to optimally-doped regime. 

PACS numbers: 74.25.Fy, 74.25. Dw, 74.72.Bk, 74.72.bk, 74.72.Dn 



The "gapless" nature of the d-wave superconductivity 
in high-Tc cuprates opened a new possibility to study the 
electronic state at low temperature through quasiparticlc 
(QP) heat transport iL2.1.A,^3.,2.A.^M.dLA2.,M- 
In the presence of disorder, QPs are created near the gap 
nodes with an "impurity bandwidth" 7 |l4j | and carry 
heat; intriguingly, those QPs are also scattered by the 
same disorder, and in some cases the effects of creation 
and scattering balance, causing the thermal conductiv- 
ity K due to QPs to become independent of the disorder 
concentration at low T. This phenomenon is called the 
"universal thermal conductivity" , which is best under- 
stood in the framework of the self-consistent T-matrix 
approximation (SCTMA) theory (iJl for BCS rf-wave su- 
perconductors. In the SCTMA theory, when 7 is in the 
universal limit k^T <^ 7 ^ Ag (Ag is the maximum 
gap), the residual component kq/T (T — s- limit of k/T) 
becomes independent of the QP scattering rate F and is 
expressed as J^] 
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with n the number of Cu02 planes per unit cell, c the c- 
axis lattice constant, and vp (^2) the QP velocity normal 
(tangential) to the Fermi surface at the node. This is 
useful because, if valid, one can obtain Ag by the bulk 
measurement of k 

Experimentally, kq/T was reported to be approx- 
imately independent of the impurity scattering in 
nearly optimally-doped YBaTCuaO^ (YBCO) 0] and 
Bi2Sr2CaCu208+5 (Bi2212) 0, and its magnitude for 
Bi2212 was consistent with Eq. Q using the angle- 
resolved photoemission spectroscopy (ARPES) results 
[l7|. Hence, it has been considered that the universal 
thermal conductivity is essentially realized in cuprates. 
However, the experimental results so far Q,y,y,|j| carry 
rather large error bars and thus it has not been clear 
to what extent kq/T is "universal". Moreover, Hussey 
et al. found an absence of kq/T in an underdoped 
YBa2Cu408, which questions the universal scenario and 



suggests a QP localization phenomenon that does not 
show up in the SCTMA theory. In this regard, an un- 
usual QP localization in magnetic field was recently ob- 
served in underdoped La2-a;Sra;Cu04 (LSCO) 0, 0| 
and, while it is obviously related to the magnetic-field- 
induced spin/charge order 1 12',' 18], the mechanism of this 
localization is not yet understood. 

Theoretically, whereas the SCTMA theory is a per- 
turbation theory, there are other non-perturbative ap- 
proaches 0, 0, 123 that tend to predict QP localiza- 
tions. Actually, when the superconducting (SC) state is 
inhomogeneous, the SCTMA theory is clearly inappro- 
priate and nonperturbative theories becomes necessary. 
Since the cuprate superconductors have an inherent ten- 
dency to develop electronic inhomogeneities in the SC 
state because of the short coherence length, and indeed 
various sorts of electronic inhomogeneities are being dis- 
covered in cuprates 0, 113, |U it is important to 
sort out the consequence of the inhomogeneity in the low- 
energy physics probed by the QP heat transport. 

In this Letter, we critically examine the universal ther- 
mal conductivity by studying YBCO, LSCO, and Bi2212 
for various doping levels; what is crucial to this work 
is that our samples are among the best single crystals 
available for these three systems, well characterized by 
various transport properties 2& i2Qj 27, 28„ ^J, and 
the K measurements are done with a small absolute un- 
certainty of less than 10%. In our experiments, we find 
three features that all indicate a breakdown of the univer- 
sal thermal conductivity upon a closer look: i) a slight 
Zn substitution for Cu induces notable suppression of 
Kg/T in underdoped and optimally-doped samples of all 
the three systems; ii) in Bi2212, both the magnitude and 
the doping dependence of the gap parameter vp/v2 ob- 
tained from Kg /T via Eq. differ significantly from the 
ARPES results jHISll; iii) a large in-plane anisotropy 
of ^2 is observed for kq/T in Bi2212. We discuss that 
the electronic inhomogeneity plays a major role in the 
departure from the universality by causing a partial QP 
localization as suggested by the nonperturbative theories. 
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FIG. 1: (color online) Zn-substitution effect on the low-T 
thermal conductivity of YBCO, LSCO, and Bi2212 in the 
underdoped and optimally-doped regimes. The doping levels 
of Bi2212 samples are indicated by the Tc values of pure ones 
(Zn-substituted ones are annealed to achieve the same oxygen 
contents as the pure ones). The solid lines are linear fits to the 
lowest-T data to extract kq/T, except for the 2 = sample 
in (f) (see text). 



High-quality single crystals of YBa2(Cui_2Zn^)30y, 
La2_a;Sra;Cui_zZn304 and Bi oSro CafCui -^Zn^)208-i-A 
are grown by a flux method |25| . a traveling-solvent 
floa ting -zone method |2^ and a floating-zone method 
HM Esl) respectively. The actual Zn concentration z is 
determined by the inductively-coupled plasma atomic- 
emission spectroscopy (ICP-AES). All the YBCO sam- 
ples are perfectly detwinned and n is measured along 
the a axis to avoid the additional electronic heat trans- 
port coming from the Cu-0 chains. For a doping- 
dependence study of Bi2212, Bi2Sr2CaCu208+5 crystals 
as well as Bi2Sr2Cai_yDyyCu208+5 (Dy-Bi2212) crystals 
(also produced by the floating-zone method) are care- 
fully annealed at 400~800°C in suitable atmospheres to 
tune the oxygen content. We label the Bi2212 samples as 
UD70, OP94, OD70, etc., by their doping regimes and 
values, and the Dy-Bi2212 samples as Dy81, Dy45, etc., 
by their values (all in underdoped regime) |32| . The 
thermal conductivity measurement in millikelvin region 
is done by a conventional steady-state "one heater, two 
thermometer" technique in a dilution refrigerator |y, lll| . 
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FIG. 2: (color online) Zn-substitution effect on the low-T 
thermal conductivity of LSCO and Bi2212 in the overdoped 
region. The data that are apparently affected by the electron- 
phonon decoupling (see text) are shown by open symbols. The 
solid lines are linear fits to estimate ko/T. 



We emphasize that our k data are reproducible with typi- 
cal variations of less than 10%, thanks to the consistently 
high quality of the crystals and a good control of the dop- 
ing level j25„ 2a, i2Sy . 

First, we show the effect of Zn-substitution on the 
low-T QP heat transport for different doping regimes. 
Figure 1 shows the n/T vs T^ plots for YBCO, LSCO, 
and Bi2212 systems in the underdoped and optimally- 
doped regimes. The T=0 intercepts of the linear fits to 
the lowest-T data give the residual QP component ko /T 
BIllllllllSIlSlIIl- Here, the striking feature is 
a strong decrease of ko/T upon slight Zn substitution 
in all these systems |33|, whereas within the SCTMA 
theory n^/T is predicted to increase with F when the 
system is not in the universal limit [isl Ib^ . We note 
that in YBCO at y = 7.00, kq/T decreases by ~20% 
with 0.6% of Zn [Fig. 1(b)]; while this behavior sug- 
gests a departure from the universal behavior given our 
small error bar (~10%) Is^l, it is not inconsistent with 
the previous results by Taillefer et al. for y = 6.9 0,0] , 
where the error bar was as large as ~40%. A more pro- 
nounced suppression in kq/T of ~40% is observed upon 
0.6% of Zn-substitution in underdoped YBCO &i y — 
6.50 [Fig. 1(a)], for which there has been no previous 
report. Apparently, the present higher-precision mea- 
surements indicate that all three systems do not strictly 
display the universal behavior in the underdoped and 
optimally-doped regimes; furthermore, the observed sup- 
pression of Hio/T by Zn substitution is opposite to what 
is expected from the SCTMA theory. If we remember 
that a non-perturbative theory that takes into account 
the gap inhomogeneities induced by impurities predicts 
a QP localization T^, it seems most likely that the ob- 
served behavior is essentially due to the Zn-induced elec- 
tronic inhomogeneity in the form of nonsuperconducting 
droplets |2^ . 

Figure 2 shows the Zn-substitution effect on k in over- 
doped LSCO and Bi2212. In overdoped cuprates, serious 
electron-phonon decoupling occurs at very low T l37l | , 
causing the QP thermal conductivity to be undetectable 
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FIG. 3: (color online) (a),(b) a-axis thermal conductivity of 
Bi2212 and Dy-Bi2212 at low T; Dy80 data are shifted up 
by 0.01 W/K^m for clarity. The solid lines are linear fits to 
extract kq/T. (c) Hole-doping dependence oi kq/T in Bi2212. 
(d) The gap parameter Vf/v2, calculated from kq/T using Eq. 
Q, as a function of Tc in the underdoped region; also shown is 
vf/v2 obtained from ARPES pl| for comparison. The dashed 
and dotted lines are guides to the eyes. 

by the conventional steady-state technique. As a result, a 
strong downturn shows up in the lowest-T data of k/T for 
overdoped samples [s^ [s^l , which is also the case with 
our data in Fig. 2. [Our optimally- doped Bi2212 also 
shows this problem below 100 mK, see Fig. 1(f)]. It is 
therefore impossible to precisely determine kq/T oi these 
samples from the present data. Nevertheless, for a qual- 
itative evaluation of the Zn-substitution effect, one can 
crudely estimate kq /T by a linear fitting to a higher tem- 
perature range, as shown in Fig. 2, and it would be safe 
to conclude that upon Zn substitution kq/T is enhanced 
in these overdoped samples. This trend is opposite to 
that observed in underdoped and optimally-doped sam- 
ples in Fig 1, and is probably understandable within the 
SCTMA theory |l5ll34l |. Hence, there is a crossover near 
optimum doping beyond which the nodal QPs behave 
more ordinarily. 

We next show the doping dependence of kq/T for 
Bi2212, where the gap parameters have been reliably 
obtained by ARPES studies [HIIJl. Figures 3(a) and 
3(b) show the low-T thermal conductivity (along the a 
axis) of Bi2212 and Dy-Bi2212 single crystals, and Fig. 
3(c) summarizes the doping dependence of kq/T. Here, 
the hole concentration per Cu, p, is determined by the 
method employing the Hall coefficient proposed in Ref. 
[38|. Clearly, kq/T decreases with decreasing doping and 
approaches zero upon entering the nonsuperconducting 
region — this behavior is essentially the same as that 
observed in LSCO [1 and YBCO The gap pa- 

rameter vf/v2 calculated using Eq. ^ is shown in Fig. 



FIG. 4: (color online) (a),(b) T dependences of a- and &-axis 
thermal conductivities for optimally-doped and underdoped 
Bi2212. (c),(d) k/T vs plots of the low-T data; the solid 
lines are linear fits to extract hq/T. 



3(d); also included in this panel are the vf/v2 values 
directly obtained from ARPES for several underdoped 
Bi2212 crystals jsJl- Not only the opposite trend in the 
doping dependence, but also the difference in the vplv2 
values from kq/T and from ARPES, demonstrate an in- 
adequacy of Eq. (Q). To be more precise, kq/T of 0.051 
W/K^m for the optimally-doped Bi2212 yields vf/v2 = 
65 and Aq = 10 meV (assuming a simple d^^-y^ ga-p); in 
contrast to vp /y?. = 20 and Aq = 35 meV obtained from 
ARPES [HBj. The reason why kq/T is much larger 
than what the ARPES-derived vf/v2 would predict via 
Eq. ^ is not clear, but this should be understood in re- 
lation with the electronic inhomogeneity, given that there 
is likely a lar ge v ariation in the SC gap magnitude in this 
material |2l|,|2J|. It might be that, depending on the na- 
ture of the inhomogeneity, sometimes the QP creation 
effect overwhelms the QP scattering effect (which may 
be the case near optimum doping) and sometimes the 
opposite happens (possibly in the underdoped regime). 
Clearly, more work is called for on the theoretical part. 

Figure 4 shows another novel property observed in 
Bi2212; namely, the QPs show very different abilities to 
conduct heat along the a and h axes that correspond to 
the two orthogonal nodal directions. Indeed, the data 
in Fig. 4 demonstrate that in the T ^ limit the 
anisotropy in kq/T is more than a factor of 2. If the 
universal scenario holds here and Eq. is valid, the 
in-plane anisotropy would mean that the SC gap is not 
the standard d^2_y2, a possibility discussed in Ref. |29l| : 
however, given the large anisotropy for T — > and the 
problems in the universal scenario, it seems more natu- 
ral to consider that the source of the anisotropy is not 
an asymmetry in the SC gap but is the breakdown of the 
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universal scenario. To corroborate this interpretation, 
the kq/T values for both a and b axes measured here 
in our optimally- doped Bi2212 are much larger than the 
values reported previously 0, ; such a dependence of 
K,o/T on the sample source suggests that the details of 
the electronic inhomogeneities matter in Bi2212. 

As for the possible source of the in-plane anisotropy in 
Bi2212, it is useful to notice that the latest "gap map" 
data obtained by the scanning tunneling microscope ex- 
periments by McElroy et al. |2ll | present clear anisot ropy 
in the electronic inhomogeneity; in Fig. 2B of Ref. |21|. 
the variation of the gap magnitude is rather smooth along 
the a axis, while it is much steeper along the h axis. If 
a stronger electronic inhomogeneity causes a larger ten- 
dency for QP localization, as the calculations based on 
the Bogoliubov-de Gennes equation suggest 0, 0| , one 
would expect Kh to become smaller than in the situ- 
ation as indicated in Rcf. 21] . Hence, the unusual in- 
plane anisotropy of kq/T in Bi2212 is likely to be related 
to the anisotropic electronic inhomogeneity, that is indi- 
rectly caused by the supermodulation structure peculiar 
to this system. 

A clear message one can take is that the QP thermal 
conductivity is unexpectedly sensitive to the electronic 
inhomogeneity, whatever the nature of the inhomogene- 
ity is; it can be the local suppression of the SC gap by 
Zn impurities [i^l, or it can be the gap inhomogeneity 
produced by dopant atoms j21| . The qualitative change 
in the Zn-substitution dependence near optimum doping 
suggests that a strong tendency to form an electronic in- 
homogeneity is present only in the under- to optimally- 
doped regime; in this regard, it is useful to note that 
the dynamics of photo-induced QPs in our Bi2212 crys- 
tals shows a sharp change at optimum doping, and a 
more BCS-like behavior is found in the overdoped regime 
[s^ . In any case, for a better understanding of the QP 
properties in cuprates, one should definitely take into 
account the QP localization due to electronic inhomo- 
geneities, and the present work demonstrates a crucial 
role of the nanoscale inhomogeneity in the low-energy 
physics of high- Tc cuprates. 
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